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Abstract 

We consider the calculation of electromagnetic fields generated by an electron bunch 
passing through a vacuum chamber structure that, in general, consists of an entry pipe, fol¬ 
lowed by some kind of transition or cavity, and ending in an exit pipe. We limit our study 
to structures having rectangular cross-section, where the height can vary as function of lon¬ 
gitudinal coordinate but the width and side walls remain fixed. For such sfrucfures, we 
derive a Fourier represenfafion of fhe wake pofenfials fhrough one-dimensional functions. 
A new numerical approach for calculafing fhe wakes in such sfrucfures is proposed and im- 
plemenfed in fhe compufer code ECHO(2D). The compufafion resource requiremenfs for 
fhis approach are moderafe and comparable fo fhose for finding fhe wakes in 2D rofafion- 
ally symmefric sfrucfures. Numerical examples obfained wifh fhe new numerical code are 
presenfed. 

PACS numbers: 41.60.-m, 29.27.Bd, 02.60.Cb, 02.70.Bf 


Submilled lo Physical Review Special Topics-Acceleralors and Beams 
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I. INTRODUCTION 


The interaction of charged particle beams and the vacuum chamber environment 
can be quantified using the concept of impedance or wakefield [ 1 ]. In order to cal¬ 
culate electromagnetic fields in accelerators several different numerical approaches 
have been suggested and implemented in computer codes. Such calculations for 
complicated three dimensional structures, however, remain a challenge even for to¬ 
day’s parallel computers. The approaches suggested in papers [2]-[4] allow one 
to obtain reliable results for rotationally symmetric and general three-dimensional 
structures even when using a personal computer. However, fully three dimensional 
calculations can be cumbersome and require long computation times. In contrast, 
for rotationally symmetric structures, the modelling is much simpler, the execu¬ 
tion time is much shorter, and the required computational resources are relatively 
modest. 

Many three-dimensional vacuum chamber components used in particle acceler¬ 
ators can be well-approximated with structures of rectangular cross-section whose 
height can vary as function of longitudinal coordinate but whose width and side 
walls remain fixed. Three examples of such geometry are shown in Fig. 1: (a) a 
corrugated structure that can be used as a dechirper [5], (b) a rectangular tapered 
collimator, and (c) a vacuum chamber in a bunch compressor. Calculation of the 
Wakefields for such structures can be greatly simplified by expanding the field gen¬ 
erated by the beam into Fourier series. As it turns out, each Fourier harmonic can be 
solved separately, and the wakefield calculated as a sum of their contributions. This 
reduces the original 3D problem into a number of 2D ones, each of which requires 
much less computing resources than the original problem. 

In this paper we derive the equations for the Fourier harmonics and propose a 
computational scheme for their solution. The scheme was implemented in a com¬ 
puter code that we call ECHO(2D). 
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FIG. 1: Structures of constant width: (a) dechirper, (b) tapered collimator, (c) vacuum 
chamber in bunch compressor. 

The paper is organized as follows. In Seetion II, we apply the Fourier transform 
in the horizontal direetion (the direetion of eonstant width), and derive a system 
of equations for the Fourier harmonies. In Seetion III, we prove that the longitu¬ 
dinal wake funetion satisfies the Laplace equation with respect to the coordinates 
of the source particle. In Section IV, we show that, in the general case, the wake- 
field of each harmonic can be described by four scalar functions of one variable. 
The description is further simplified in Section V where we assume that the entire 
structure has a horizontal symmetry plane, and then analyze the longitudinal and 
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transverse Wakefields near the system axis. An important symmetry relation for the 
wake funetion is established in Seetion VI. In Seetions VII and VIII we present a 
numerieal algorithm for solving equations derived in Seetion II. Several numerieal 
examples obtained with the eomputer eode ECHO(2D) are deseribed in Seetion IX. 
The results of the paper are summarized in Seetion X. 

II. PROBLEM FORMULATION AND FOURIER EXPANSION 



FIG. 2: Charged particle bunch moving through an accelerating structure. 

We eonsider a buneh of partieles moving with the veloeity of light c through 
a metallie strueture eomprising an ineoming pipe, a transition or eavity, and an 
outgoing pipe, an example of whieh is shown in Fig. 2. The buneh is eharaeterized 
by the eharge distribution p and the eleetrie eurrent density j = cp. We assume 
that the buneh is moving along a straight line parallel to the longitudinal axis of 
the system, and we negleet the influenee of the wakefields on the buneh’s motion. 
Our problem is to find the eleetrie and magnetie fields E, H, in the domain fi 
whieh is bounded transversally by perfeetly eondueting metallie walls, defined by 
the boundary dQ of the domain Q. We need to solve Maxwell’s equations with the 
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boundary conditions: 


VxH = —D + j, 
dt ■' 

V-D=p, 

B=pH, 


E(t - 0) - Eq, 


nx E = 0, X G dQ, 


VxE = -^B, 
dt 

V 8 = 0, 

D = eE, r e Q, 

H{t = 0) = Ho, re Cl, 

( 1 ) 


where ^o, Hq, is an initial eleetromagnetie field in the domain Q and « is a unit 
vector normal to the surfaee dQ.. While we eonsider the propagation of beams 
through vaeuum, for the generality of the eomputational algorithm we ineluded in 
Eqs. (1) arbitrary dieleetrie eonstant e and magnetie permittivity p. 

We ehoose a coordinate system with y in the vertical and x in the horizontal 
direetions; the z eoordinate is direeted along the longitudinal axis of the system. 
The structures eonsidered in this paper have constant width 2w and a fixed loeation 
of the side walls (in the x-direetion). 

The charge density is localized within the interval 0 < x < 2w and vanishes at 
the ends of the interval, p(0, y, z) = p(2w, y, z) = 0. Such a function can be expanded 
in Fourier series 


^ oo 

p{x,y,z) = - ) Pm{y,z)sm{kx^mX), 
w 


n 


Pmiy,z) 




m=l 

•2w 


kx,m ~ ^ ^5 

2w 


p{x, y, z) sm(kx,mX)dx. 


( 2 ) 


It follows from the linearity of Maxwell’s equations, Eqs. (1) and Eq. (2), that the 
eomponents of eleetromagnetie field can be represented by infinite sums: 


COSi^kx^mX)- (3) 


Hx 


Hx,m 


( \ 
Ex 


( \ 
p 

^x,m 

Ey 

II 

Ms 

^y,m 

Cm{kx,mX), 

Hy 

II 

Ms 

^y,m 


m=l 

y^z,m ) 


V_ 

m=l 
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For each mode number m we can write an independent system of equations 



(4) 


Hence we have reduced the 3D problem, Eqs. (1), to an infinite set of independent 
2D problems, Eqs. (4). Our approach is similar to one used in [6] for solving 
Maxwell’s equations in the frequency domain. 

III. LONGITUDINAL WAKE AS A HARMONIC FUNCTION 

Let us consider a line-charge beam with vanishing transverse dimensions, 
p(xo, yo, X, y, s) = Q6{x - xo)6(y - yoU(s), 


jz(xo, yo, X, y, s) = cp(xo, yo, x, y, s). 


(5) 


where xq, yo, define the transverse offset of the beam, s = z - ct is the local longi¬ 
tudinal coordinate in the bunch, Q is the bunch charge and /l(j') is the longitudinal 
bunch profile [for a point charge, /l(5) = d(5)]. The longitudinal wake potential W\\ 
at point {x, y, s) is defined as [ 1 ] 



( 6 ) 
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where the eleetrie field on the right-hand side is the solution to Maxwell’s equation 
with the sourees of Eqs. (5) (this field, of eourse, is also a funetion of Xq and yo 
omitted in the arguments of for brevity). 

It is well known [7] that the longitudinal wake potential satisfies the Laplaee 
equation for the eoordinates x, y, of the witness partiele 



(7) 


In order to express the wake potential in struetures of eonstant width through 
one-dimensional funetions we need to prove that the longitudinal wake potential 
satisfies also the Laplaee equation with respeet to the eoordinates xq, yo, of the 
souree partiele. This statement will aetually be proven below for arbitrary strue- 
tures without any restrietions on the strueture geometry. The only requirement im¬ 
posed on the system is that the ineoming and outgoing waveguides have perfeetly 
eondueting walls. 

The proof is based on the direetional symmetry relation between the wake po¬ 
tential Vk|| for the ease when the partiele travels through the system in the positive z 
direetion, and the wake potential eorresponding to the ease when the partiele 
travels through the same system but in the negative z direetion. It was shown in 
Ref. [8] that the two wake potentials are related by 


IT|] \r 2 , ri,s)- Wiiiri,r 2 , s) = 2 cZe{rur 2 )X{s), 


( 8 ) 


where 



(9) 


eo is the permittivity of vaeuum, r\ = (xq, yo) and r 2 = (x, y) are two dimensional 
veetors, and (pA, 4 >b, are the Green funetions for the Laplaeian in the ineoming and 
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outgoing pipes with cross-sections, Sa,Sb- 


A 0 ^(r;, r) = -Co - r,), r e Sa, 

(f>A(ri,r) = 0, reds A, 

A(f>B(ri, r) = -Co - r,), r e Sb, 

(pB{ri,r) = Q, redSB, i = 1,2, 


( 10 ) 


where dS is the boundary of S. Note that r\ and r 2 are offsets of, respectively, the 
leading and the trailing particles in W^|, while r 2 is the offset of the leading particle, 
and r\ is the offset of the trailing particle in Wy”’- Note also that, according to the 
general property of the wake potentials (7), satisfies the Laplace equation 

I d^ d^ \ (_) 

^ r^ll ^0’^^0, ^) = 0. (11) 

\^jcq 

With the help of Green’s first identity (n is the outward pointing unit vector 
normal to the line element dl). 


V(p ■ Vtf/ds 

we can rewrite Eq. (8) as 


J V(p ■ Vipds = I ifdnij/dl - I (pAiffds, 
s Jds Js 


\r2,ru s) - Wi\iri,r2, s) = 2A(s)[(f>Birur2) - (f>Airur2)], ( 12 ) 


where we have used the symmetry of Green’s functions, 0 a(^ 1 ,^ 2 ) = 0a(^ 2 ,^ 1 ), 
0b(^i, ^ 2 ) = 0b(^ 2, ^i)- We now apply the Laplacian operator to the coordinates Xq, 
yo in Eq. (12) and use relation (11) to obtain 


d^ d^ 

77:2 + ^ I y ^ 


dxi 




(13) 


This proves that the wake potential, in addition to being a harmonic function with 
respect to the coordinates of the trailing particle, is also a harmonic function with 
respect to the coordinates of the leading particle. To our knowledge, this general 
property of wakefields has not been previously reported in the literature. 
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IV. WAKE POTENTIAL REPRESENTATION IN STRUCTURES OE CONSTANT 
WIDTH 


It is well known [7] that for rotationally symmetric structures the wake potential 
can be represented through one dimensional functions, with only one function for 
each azimuthal mode number m, 

oo 

lV|l(;co, i/o, X, y,s)='Y^ W^isy^r"^ cos{m{e - Oo)), 

m=0 

xo = ro cos(0o), i/o = ro sin(0o), x = rcos(6), y = rsm(6). (14) 

As it turns out, a similar, comparably simple representation exists for the wake 
potential in structures of rectangular cross-section and constant width; however, in 
the latter case four one-dimensional functions are needed for each mode number m. 
Substituting Eq. (5) into (2) we find for the charge distribution 

1 

p(xo, yo, X, y,s) = - 
w 

Prniyo, y, s) = Q6{y - yo)^is). (15) 


'^Pmiyo, y, s) sinik^^^xo) 


It then follows from Maxwell’s equations (4) that all components of the fields will 
be proportional to sin(k;c,m-^o) and using Eqs. (3) and (6) we find 


^ oo 

W||(xo, yo, X, y,s) = - V lV,„(i/o, y, s) sm{k:c,mXo) sin(k_,_„,x). 

w ^ 

m=\ 


(16) 


If we insert representation (16) into Eq. (7), we obtain the one-dimesional Helmholz 
equation 


dy 


-^Wrniyo, y, s) - (kj,^m) Wmiyo, y, s) = 0, 


which has a general solution of the form 


(17) 


Wn,(yo, y, s) = W^„(yo, s) cosh{kjc,my) + (?/o, s) sinh(k^,,„j/). (18) 
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Finally from Eq. (13) it follows that functions W^^{yo, s) and W^„(yo, s) satisfy equa¬ 
tions 


52 


9yo 


^W:,(yo,s)-(h^J^W^Jyo,s) = 0, 




2 " m 


dy\ 

These equations ean again be easily integrated with the result: 

^ 00 

lT|l(;co, yo, = ~'Yu y^ sin(fc;,,„,A:o) 

where 


(19) 


( 20 ) 


m=l 


M^m(?/o, y, S) = [Wl"{s) eosh(fc;,,^i/o) -H W'^{s) sinh(fc;,,„i/o)] co&KKmy) 

+ [1T,);"(5) co^hik^^myo) -H Wl'{s) sinh(fc^,^^o)] smh{k^^,ny). (21) 


Thus, we have proven that, in struetures of eonstant width, for eaeh mode number 
m four funetions are needed to eompletely deseribe the longitudinal wake potential. 
These funetions ean be ealeulated as follows 


1T^‘^ = 1^^(0,0, ^), 
1 d 


1 d 

= ^^W„{0,0,s), 


’ ' w 


1 T„,( 0 , 0 , 5), 


IT" 

' ' vy. 


kx ffi Ox 

1 


kx,m “ ikx,m)^ dxdxo 

where the m* modal eomponent of the wake potential 

X oo 

[E,^,„{y,z, t)]t=(,-x)/cdz 

00 

is exeited by a eharge distribution that does not depend on x. 


W„,{0,0,s), (22) 


(23) 


Prniyo, y, s) = Q6{y - yoUis). 


(24) 


With a knowledge of the longitudinal wake we ean ealeulate the transverse 
wakes. For example, the vertieal wake potential, Wy, ean be easily found through 
the Panofsky-Wenzel theorem [7] 


d d 

— Wy(xo,yo,x,y,s) = —Wii(xo,yo,x,y,s). 


(25) 
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In the next seetion we will analyze both the longitudinal and the transverse wakes 
assuming that the strueture under eonsideration has a symmetry axis. 


V. STRUCTURES WITH HORIZONTAL SYMMETRY PLANE 

Let us eonsider a strueture of eonstant width 2w that also has a vertieal symmetry 
plane, aty = 0. Struetures in Fig. 1 (a) and (b) possess this symmetry; henee, they 
have a symmetry axis loeated at x = w, y = 0. Due to the symmetry, the wake 
potential satisfies the equation 

W^l(xo, yo, X, y, s) = W||(xo, -yo, x, -y, s), (26) 

and Eq. (21) simplifies: 

W,n{yo,y, s) = W^^(s)cosh{kjc,myo)cosh(k^^,„y) + W'J(s)smh(k^^,nyo)smh(k^^^y). 

(27) 

Note that 

Wmiyo, y, s) = W^iy, yo, s). (28) 

Let us eonsider the transverse wakes in sueh struetures. We first introduee the 
integrated wake funetions (sometimes ealled the step funetion response) 

s‘: = f s:: = f w-(/)j.'. (29) 

kJ—OO kJ—OO 

It then follows from (25) that the transverse wake funetion ean be written as 

1 “ 

Wy{xo, yo, X, y,s) = - V k^^„Wy^m{yo, y, s) sin(k_,_„,xo) sm(kj,^^x), (30) 

w 

m=\ 

where 

Wy,miyQ, y, S) = cosh(k^^^yo) smh(k^^^y) + S sinh(k;,,,„i/o) cosh(k^^^y), 

Wj,(xo, yo, X, y,s) = - V k^^^W^^miyo, y, s) ?,m{k^^rnXo) eos(kx,mJc), (31) 

w ^ 

m=l 

W^^miyo, y,s) = Scosh(k^^^yo) cosh(k^^^y) + S smh(k^^^yo) smh(k^^rny)- 
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Representations (30), (31), are valid for arbitrary offsets of leading and trailing 
partieles. 

For small offsets near the symmetry axis, x = w, y = 0, the transverse wake 
potential is usually expanded in Taylor series. 


(9 d 

Wy(w, yo, w, y, s) « yo-^Wy(w, yo, w, 0, + y— Wy(w, 0, w, y, 5)1^^^^. (32) 

The first term in (32) is usually ealled the transverse dipole wake in the i/-direetion. 
It ean be ealeulated as follows 


d 1 ^ 

Wy,As) = ^Wy{w,y^,w,Q,s)\ „ = - V (33) 

oyo w 

m=l ,odd 

The seeond term in (32) is ealled the transverse quadrupole wake in i/-direetion; it 
is obtained by 

d 1 ^ 

" m=l,odd 


The transverse wakes in the x direetion are obtained by equations eorresponding to 
those of Eqs. (33), (34). Note that Wy^q{s) = -W^^q^s). 

In numerieal ealeulations of struetures with symmetry we ean use the approaeh 
of paper [9] that allows us to reduee the ealculation domain in half. Indeed the 
eharge distribution (24) ean be written as a sum of symmetrie and antisymmetrie 
parts 


pmiyo, y, s) = pliyo, y, s) + p^iyo, y, s). 


(35) 


where 

Pmiyo, y,s) = ^Q [6(,y - yo) + 5{y + yo)] A(s), (36) 

Pmido, y,s) = [6{y - yo) - 5{y + yo)] A(s). (37) 

In problems with the symmetrie driving eharges (36), the tangential eomponent of 
the magnetie field will be zero in the symmetry plane (the so ealled “magnetie” 
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boundary condition). In problems with the antisymmetric driving charges (37) the 
tangential eomponent of the eleetrie field will be zero in the symmetry plane (the 
“eleetrie” boundary eondition). Thus, instead of solving the system of equations 
(4) in the whole domain, one ean solve two independent problems in half of the 
domain: one problem with the “magnetie” boundary eondition at ^ = 0 and one 
problem with the “eleetrie” boundary eondition at y = 0. This is true not only 
for the line-eharge eurrent distribution (5), but for any arbitrary three dimensional 
eharge distribution p(v, y, z, t). From solutions Wj^iyo, y, s) and (i/o, y, s) of the 
two problems we ean easily find the one dimensional modal funetions in Eq. (27): 




(38) 


VI. SYMMETRY RELATIONS FOR THE TRANSVERSE WAKE POTENTIAL 

In the general ease, from the direetional symmetry relation, Eqs. (8), diseussed 
in Seetion III, one eannot immediately infer the eorresponding symmetry for the 
transverse wakes taken in the positive and negative direetions. It turns out, how¬ 
ever, that for struetures of eonstant width with a symmetry axis, sueh a direetional 
symmetry relation for the transverse wake potential ean be proven. Indeed, from 
the symmetry relation (28) and Eq. (20) follows the symmetry of the longitudinal 
wake potential. 


(39) 


W||(xo, yo, X, y, s) = Wi\{x, y, xo, yo, ^)- 


Henee the direetional symmetry relation (12) for the longitudinal wake potential 
ean now be written as 


Wj \rur2,s)- Wi\{ri,r2,s) = 2T(5) [^^(ri, r2) - 0A(ri, ^ 2 )]; (40) 
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and from the Panofsky-Wenzel theorem (25) it follows that 


W[ \ri,r 2 , s) - Wx(ri, r 2 , s) = 2A(s) [V^sCrj, r 2 ) - V^aC^i, ^ 2 )], 


A(^) = 



A{s')ds'. 


( 41 ) 


In the last equation we have used the veetorial notation for the transverse wake, 
Wx = (W„ W,). 

Let us now ealeulate the direetional relation for the dipole and quadrupole wakes 
defined by Eqs. (33), (34). Denote the vertieal aperture of the ineoming reetangular 
pipe by Ig^ and that of the outgoing pipe Igs- (Both pipes have the same width 2w.) 
The Green funetion for Poisson’s equation in a reetangular pipe, Do, y), has 
been obtained by Gluekstern, et a/ [10]. For xq = w, their result reads:^ 


(pAiw, yo, x,y) = - 

neo 


^ 00 

--y 

TT(=rs 


ST 


e cosh 


2g 


n=l 


neosh^ 


. nn nn^ 

sin —(y + gx) sm —(yo + gA) 
2g 2g 


+ \n\\{x - wf + {y - yof]- 

I L Jo- 


sinh 


2 n{x-w) 


— + eos^ 

Aqa 49a 


{y + yoY 


sinh 


2 n{x-w) 


49 . 


In [(x - wf + {y- yof\. 


(42) 


The last term is singular in the limit x ^ w and y ^ yo- All other terms are 
finite for all admissible x and y, and the sum in the first term rapidly eonverges. 
The Green’s funetion of the outgoing pipe, ^b(«;, yo, x, y), is obtained from (42) 
by replaeing gA gs- It is then easy to see that the singular terms eaneel in the 
differenee (pA-(pB\ henee this funetion is finite within the domain of its definition. If 
we follow the approaeh of paper [11] and use Eqs. (39), (40), (45) from that paper. 


^ There is a typo in Eq. (5.11) of Ref. [ 1 0] that is corrected in (42). 
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then from Eq. (41) our final result reads 


[^yM - 

J(^) - Wy,gis) = 2cA(s) [Zy^yigA) - Z,,,(^fi)]. 
nZo 


Zy^iCj) — 
Zy^qig) — 


12^2 

ttZq 

24^2 


1 + 




m 


m=l 




m=l 


^ g2nm{w/g) 

2m - 1 


^ g7T{2m-l){w/g) 


(43) 


VII. TE/TM SCHEME IN MATRIX NOTATION 


In this seetion we deseribe a partieular realization of the implieit TE/TM seheme 
introdueed in [4]. The seheme will be diseussed in the eontext of the Einite Integra¬ 
tion Teehnique [12]. 



FIG. 3: One cell of cell complex G showing the positions of the voltage and magnetic flux 
components 

We will eonsider Maxwell’s equations in their integral form on a domain Q with 
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the linear non-dispersive eonstitutive relations: 


(£Edl = (ffi Bds, (f 

) Hdl = ^ ^ Dds + ^ jds. 

V5 c Q, 

h s k 

s s 



= (|y|) pdv, Bds = 0, 

vy c a 

dV 

V dV 



D = eE, B = iiH, 

Vx 6 kl. (44) 


Let us start by introdueing a grid-based deeomposition of the entire eomputa- 
tion domain Q into eell eomplex G. We use here a three dimensional Cartesian 
mesh in x, y, z eoordinates, with the eorresponding indexing i, j, k. Unlike in finite- 
differenee methods where one starts by alloeating the field eomponents, we begin 
by taking the voltage along eell edges and the magnetie flux through eell faeets as 
the eomputational unknowns: 




I Edl, bff = II Bds, 

L§- \}\J 


(45) 


where ?? is a mesh multi-index and L^,S& 6 G. Solving Faraday’s law in integral 
form for the front surfaee shown in Fig. 3 yields: 

d 


“ 1 “ ^y,ij,k ^x,i,j+l,k ^y,i+\J,k ~ 


dt 


-byi 


z,U,k- 


Note, that this representation is still exaet, as is (by definition) the exaet eleetrie 
voltage along one edge of the eell, and similarly b§^ represents the exaet value of the 
magnetie flux density integral over the eell surfaee. If we eompose eolumn veetors 


(46) 


from all voltage and flux eomponents, we ean write the eombination of all equations 
over all surfaees in an elegant matrix form as 

d. 



( \ 

ex 



e = 


b = 

by 




bz. 


Ce = 


dt 


b. 


(47) 


17 







The matrix C picks the affected components out of the long vector to make up the 
corresponding equation. C is thus the discrete curl operator over the mesh G. With 
an appropriate indexing scheme the curl matrix has a 3x3 block structure: 


C = 





(48) 


The double-banded, topological P;^, Py, P^-matrices take the role of discrete partial 
differential operators. 

The second important differential operator in Maxwell’s equations (1) is the di¬ 
vergence operator. In order to construct a discrete divergence operator we integrate 
Maxwell’s equation Bds = 0 over the entire surface of a mesh cell depicted in 

dV 

Fig. 3. By adding up the six fluxes for each cell and by writing down all such equa¬ 
tions for the entire mesh we obtain a discrete analogue to the divergence equation: 


Sb=0, S = (P, Py P,). 


(49) 


The discretization of the remaining Maxwell equations requires the introduction 
of a second cell complex G which is dual to the primary cell complex G. For the 
Cartesian grid the dual complex G is defined by taking the foci of the cells of G as 
gridpoints for the mesh cells of G. We again introduce the computational unknowns 
as integrals 

f jds, (50) 

where is a mesh multi-index and L^, e G. 

Following an equivalent procedure for the remaining Maxwell’s equations, but 
using the dual cell complex G, we obtain a set of four discrete equations represent- 


= I Hdl, d^ = 11 Dds, 


j& 
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ing Maxwell’s equations on a dual grid: 


Ce = -|b, C-h=|d.j, 

Sb = 0, Sd = q, e, b 6 G, h, d e G, 

S = (-P. -p; -p:), (51) 

where the asterisk denotes the Hermetian adjoint operator. 

Equations (51) are eompleted by the discrete form of the material relations (con¬ 
stitutive equations) which appear (in the simplest linear case) as matrix equations 


e = cMf-id, h = cM^-ib. (52) 

with Mf-i the discrete inverse permittivity matrix, and M^-i the inverse permeability 
matrix. 



FIG. 4: Grid doublet (G, G) in {y, z)-plane. 

In the case of Cartesian grids (or, more generally, whenever the primary and dual 
grids are orthogonal) all material operators can be represented by diagonal matrices. 
Note that the material matrices contain both averaged material parameters, and the 
lengths and areas of the grid edges and faces, respectively. In the standard staircase 
approximation the material matrices contain elements (without double indices for 
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FIG. 5: Curved PEC-boundary in a Cartesian mesh. 


simplicity of notation) 

^pijk ~ ^ J ^pijk ~ (^/^) ? ’ (53) 

b pijk ^pijk 

with p = x,y, z, and the face areas and edge lengths of the primary and seeondary 
grid given by 5, L, 5, L, respectively. 

In the eonformal seheme [14], we allow the eells of the computational grid to be 
only partially filled by a perfeetly eleetrie eonducting (PEC) material and with an 
arbitrarily shaped interface. To model such a case we modify only the elements of 
the material matriees: 



= (C6)-' 


^pijk 

S pijk 


fj^p]jk = W 


^pijk 

Lpijk 


(54) 


where s, I denote the reduced eell areas and lengths, ineluding only those parts 
inside the eomputational domain (outside PEC material), as shown in Eig. 5. 

The complete set of equations (51) and (52) is referred to as Maxwell’s grid 
equations. 

Eor our 2D rectangular structures. Maxwell’s equations (1) reduce to the modal 
equations (4). In addition, the dual three-dimensional cell eomplex (G, G) reduees 
to the two-dimensional plane eomplex shown in Eig. 4. In Maxwell’s grid equations 
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the discrete operator reduces to the diagonal matrix k;^ = kx^mi-, where I is the 

unit matrix. 

System (51) is a time-continuous and space-discrete approximation of problem 
(4). The next step is a discretization of the system in time. The field components can 
be split in time and the leap-frog method can be applied. With “electric/magnetic” 
splitting, a well known Yee’s scheme [13] will be obtained. However Yee’s scheme 
has large dispersion errors along the grid lines. In the following we consider an 
alternative TE/TM splitting scheme [4], one that does not have dispersion errors in 
the longitudinal direction. 

Let us rewrite Eqs. (51) in an equivalent form 


d 


dT 

where t = ct and 


u = Mu(TuU-r Lv-ju), 




V = Mv(TvV-L*u-jv), 


0 0-1 
0 0 k, 


0 0-1 
0 0 k, 
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/- - 

o 

V _ 



0 P. 0 


-P, 0 0 
0 0 0 


(55) 


(56) 


Applying TE/TM splitting [4] of the field components in time to system (55), 
the following numerical scheme is obtained 

j.n+0.5 _ j,«-0.5 / j,«+0.5 i |,n-0.5 \ 

M„ I T„-^-+ Lv" - ju” , 


At 


2 

rrH+l I 


-I- V' 

^ ’ _ -XM I'T' ’ ^ T i n+0.5 

At 


L*u«+o.5 j^ 


(57) 
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Two-layer operator-difference scheme (57) acquires the canonical form [15] 


v”+' - v" 

^ +Ay” = r, 
At 


where 


B 


Mu ' - crTu 0 

2aL* Mv~^ - aTv 


y = 


u 


«-0.5 



( 

\ 

, A = 

-T 

A u 

-L 

/ 


1 

H 

< 


£/1 _ 




1 

Ju 

;h+0.5 
\Jv J 


a = 0.5At. 


It was shown in [4] that stability of the method is insured if 

Q = B - aA > 0. 

If the matrix Q is positive definite, then we can define a discrete energy as 

lT" = 0.5<Qy",y”) 

and the discrete energy conservation law holds 


At 




«+0.5 2/7+0.5 


), 



ef + e" 


"j«+0.5' 

Jx 

g«+0.5 ^ Q_5 

«+l ^ n 

y y 

'■/I+0.5 _ 

’ J 

j«+0.5 


„«+0.5 1 p«-0.5 


i” 


( 58 ) 


(59) 


(60) 


( 61 ) 


The stability condition (59) in the staircase approximation (of the boundary in vac¬ 
uum) can be rewritten in the form 


I-(r2(Az)-2p^P* >0. 


(62) 


This last condition resembles the well-known stability condition of the explicit 
FDTD scheme for one-dimensional problems. The maximal eigenvalue of the dis¬ 
crete operator P^P* fulfills the relation A,nax < 4, and the stability condition of the 
scheme in the staircase approximation of the boundary in vacuum reads 


At < At. 


(63) 
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On an equidistant mesh the implieit seheme (57) has a seeond order loeal approxi¬ 
mation error in homogeneous regions. 

Relation (63) does not eontain information about the transverse mesh. Henee the 
transverse mesh ean be ehosen independent of stability eonsiderations. Following 
the eonventional proeedure [4], the dispersion relation ean be obtained in the form 



where O = 0.5mAT, Ky = O.SkyAy, = O.Sk^Ax. With the “magie” time step 
At = Az, the seheme does not have dispersion in the longitudinal direetion. 

Following the approaeh of [4] it is easy to show that eharge eonservation holds: 



where funetions with overbars are defined as 



VIII. TE/TM SCHEME IN COMPONENT NOTATION 

In this seetion we rewrite seheme (57) in eomponent notation and diseuss the 
algorithm of numerieal solution. 

The field eomponents ^ at time level n + 0.5 ean be found as 



( 64 ) 
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where h", h” are auxiliary veetors and operator matrix W,, is a sparse one, 


I + a" . 


( 65 ) 


The field eomponents at time level n + 1 ean be found as 







( 66 ) 


where are auxiliary veetors and operator matrix W/, is a sparse one, 


W, = I + + M^-.P,M,-.P;]. 


(67) 


In the stairease approximation of the boundary, the material matriees M^-i, M^-i 
are diagonal, and the operators (65),(67) are tri-diagonal matriees. The equations 
involving these operators ean be solved easily by the “elimination” method [15]. 
However, the stairease approximation of the boundary results in first order eonver- 
genee in the wake potential. In order to obtain seeond order eonvergenee and avoid 
time step reduetion, we use eonformal method in the same way as deseribed for the 
rotationally symmetrie ease in [16]. 

In the original problem formulation (1), we assumed that the strueture walls are 
perfeetly eondueting. However, in order to use the Fourier expansion (3), it is suf- 
fieient to require only that the walls at x = 0 and x = 2«; be perfeetly eondueting. 
The bottom and top walls of the strueture ean have different boundary eonditions; 
for example, they ean be speeified to be metallie with finite resistivity. A similar 
approaeh for the ease of rotationally symmetrie geometry and the use of the stair¬ 
ease approximation for the boundaries was deseribed in [17]. In ECHO(2D) we 
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have realized a eonformal version of the method for ineluding boundaries with fi¬ 
nite eonduetivity. Details of the implementation of our algorithm will be published 
elsewhere. 

IX. NUMERICAL EXAMPLES 

In this seetion we eonsider two example ealeulations: a eorrugated pipe 
(“deehirper”) with perfeetly eondueting walls, and a tapered eollimator with re¬ 
sistive walls. We eompare results obtained with ECHO(2D) with those of several 
3D eodes. 

For our deehirper example we used the parameters of the deehirper experiment 
performed at the Pohang Aeeelerator Laboratory [18]. The geometry is shown in 
Fig. la. The deehirper strueture is 1 m long, has a width 2w = 50 mm, and is made 
of aluminum. The eorrugations are eharaeterized by the period p = 0.5 mm, height 
h = 0.6 mm, and aspeet ratio h/t = 2. The gap between the jaws is adjustable; at a 
nominal setting the full gap 2a = 6 mm. 

We have performed ealeulations of the wakes for a Gaussian buneh with rms 
length cr^ = 0.5 mm. The transverse dimensions of the buneh were assumed negligi¬ 
bly small. In the right plots of Figs. 6-7, we show some of the harmonies W^(k^, s) 
and s), ealeulated for 2w = 5 em. If the buneh passes through the mid¬ 

dle of the strueture, at xq = w, then only odd harmonies are exeited. In the left 
plots of Figs. 6-7, these funetions are shown as eontinuous funetions of kx (denoted 
with subseript m removed); they are independent (see Appendix) of the width pa¬ 
rameter 2w. (Here we have used 2«; = 10 em to obtain a denser sampling of the 
funetions in the figures.) Finally, in Fig. 8, we eompare the longitudinal (left plot) 
and the dipole (right plot) wakes obtained with ECHO(2D) with the existing, fully 
3D eode, ECHO(3D) [4]. In the ealeulation of the wake potentials, where the width 
2w = 5 em, we used only the first 15 odd harmonies. We see that the agreement 
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with the existing eode is very good. 



V 



FIG. 6: For the Pohang dechirper: function s) (left), and the harmonics m = 1, 9, 

19, of W^is) (right). 



FIG. 7: For the Pohang dechirper: function k^W^^kx, s) (left), and the harmonics m = 1,9, 
19, of kl^W^^is) (right). 

In Table I we eompare the calculated average wakes (the loss factor in the lon¬ 
gitudinal case and the kick factors in the transverse cases) with numerical results 
obtained with the time-domain code T3P of the code suite ACE3P [19]. The results 
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FIG. 8: Longitudinal wake potential (left) and dipole wake potential (right), calculated from 
first 15 odd harmonics. 

in the table, denoted by r, have been normali z ed to analytieal results, valid for in¬ 
finitesimally small eorrugations (see [18]). We see good agreement exeept for the 
case of the dipole wake. (Note, however, that we have good reason to believe that r 
should always be < 1.) Incidentally, it should be noted that in [18] the strength of 
the wakes of the Pohang dechirper were measured, both the longitudinal and trans¬ 
verse wakes, and it was estimated that their measured (effective) r was about 0.75, 
a result that is in good accord with the ECHO(2D) results. 

TABLE I: For the Pohang dechirper, for a Gaussian bunch with = 0.5 mm: the ratio 
of the numerically obtained and the analytical result(see [18]), that we denote by r. We 
compare the results given in Ref. [19], r[T3P], with the new results, r[ECHO]. 

Wake r[T3P] r[ECHO] 

Eongitudinal, loss factor 0.84 0.83 

Dipole, kick factor 1.08 0.79 

Quad, kick factor 0.73 0.73 

As the second example problem we consider a symmetric, tapered collimator 
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FIG. 9: Longitudinal wake potential of tapered collimator. 

(see Fig. lb). The dimensions are: width 2w = 10 em, height of large pipe 
2a = 10 cm, length of tapers T = 5 cm; height and length of the minimum gap 
section, 2b = 2 cm and L = 12 cm. As was mentioned above, ECHO(2D) is capa¬ 
ble of modelling structures with metallic walls of finite conductivity. The tapered 
walls and the walls of the minimum gap section are taken to have conductivity 
cr = 100 S/m (a poor conductor), while the remaining surfaces are assumed to be 
perfectly conducting. In the left plot of Fig. 9 we compare the longitudinal wake 
for this collimator with one that has the same geometry but is perfectly conduct¬ 
ing. The Gaussian bunch in the simulations has an rms length cr^ = 0.25 cm. Both 
wakes were obtained with ECHO(2D). In the right plot of Fig. 9 we compare the 
ECHO(2D) wake potential with the one obtained using a fully three-dimensional, 
commercially available code CST [20]. The good agreement between the results in¬ 
dicates a good accuracy of the conformal meshing and the resistive wall modelling 
in ECHO(2D). 

Finally, in order to demonstrate the capabilities of the new code to calculate 
wakes of very short bunches, we give in Fig. 10 the calculated loss factor x = 
(VV]|(0,0, s)) as function of bunch length in the 1 m-long Pohang dechirper (black 
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FIG. 10: Loss factor for very short bunches in the dechirper (in black). The grey symbols 
give the difference between the loss factors for 2 m-long structure and a 1 m-long structure. 

plotting symbols). Note that the shortest bunch length used in the calculations is 
CTj, = 10 jjm. For the shorter bunches, the loss factor in a 1 m-long structure is 
larger than the steady-state (periodic) result. We can see that by plotting on the 
same figure the difference of x for a 2 m structure minus that for a 1 m structure 
(the gray curve). Note that the analytical, asymptotic value of the loss factor, for 
CTj, ^ 0, is X = ZocL/(27ra^) • (x^/16) = 1234. V/pC [we have taken Zq = 377 Q, 
structure length L = 1 m, and half-gap a = 3 mm], which agrees well with the 
linear extrapolation of the ECHO(2D) steady-state results (the gray curve). 


X. SUMMARY 

In this paper we presented a new method for solving electromagnetic problems 
and calculating Wakefields excited by a relativistic bunch in a structure that is char¬ 
acterized by a rectangular cross-section whose height can vary as function of longi¬ 
tudinal coordinate but whose width and side walls remain fixed. Using the Fourier 
expansion of the fields, currents and charge densities, we derived a Fourier repre- 
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sentation of the wake potential in terms of four one-dimensional funetions for eaeh 
harmonie. We proved that the longitudinal wake is a harmonie funetion with respeet 
of the eoordinates of the leading eharge. For a strueture that has a horizontal sym¬ 
metry plane we also proved an important symmetry relation of the wake funetion. 

A numerieal method was proposed for solving Fourier harmonies of the fields. 
The method does not generate dispersion in the longitudinal direetion and it eon- 
serves the energy and eharge in the ealeulations. The eomputer resourees required 
to solve the system with several tens of the lowest harmonies are moderate and 
eomparable to those needed for 2D rotationally symmetrie ealeulations. 
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Appendix A: Universal continuous functions for structures of constant width 

The modal eharge distribution, Eq. (24), in Maxwell’s equations, Eqs. (3), is in¬ 
dependent of the strueture width, 2w. The same holds for the boundary eonditions. 
Henee the modal solution depends on 2w only through the harmonie number 
This faet allows us to eonsider eomponents of the eleetromagnetie fields as eontin- 
uous funetions of the harmonie number, k^. The same is true for the longitudinal 
wake potential. If the eontinuous funetion W\\(yo,kx,y, s), k^ 6 [0, oo), is known, 
then we ean find the solution for any arbitrary width 2«;i as 




m=l 


(Al) 

In order to find the funetion W\\{yQ,k^,y, s), k^ 6 [0, oo), we ean proeeed in 
the following way. We take width 2w large enough in order that the sampling of 
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W\\{yQ,kx,y, s) with step kx,m+i - kx,m = n/(2w) is dense enough for smooth inter¬ 
polation of the function between the discrete values kx^m = nm/(2w). Due to the 
existence of perfectly conducting walls at ;c = 0 and x = 2w, the discrete Fourier 
harmonics do not contain a “zero” harmonic with kxfi = 0- At the same time it 
is reasonable to find the solution to this problem as well. Then, using interpola¬ 
tion, we obtain a universal continuous function W\\(yo,kx,y, s), kx 6 [0, oo), which 
is independent from the halfwidth w. 
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